Collective modes of Fermi superfluid containing vortices along the BEC-BCS crossover 



in 
o 
o 

(N 
o 

Q 

m 
(N 



Tarun Kanti Ghosh and Kazushige Machida 
Department of Physics, Okayama University, Okayama 700-8530, Japan 
(Dated: February 2, 2008) 

Using the coarse-grain averaged hydrodynamic approach, we calculate all low energy transverse 
excitation spectrum of a rotating Fermi superfluid containing vortex lattices for all regimes along 
the BEC-BCS crossover. In the fast rotating regime, the molecular BEC enters into the lowest 
Landau level, but the superfluid in the unitarity and the BCS regimes occupies many low-lying 
Landau levels. The difference between the breathing mode frequencies at the BEC and unitarity 
limit shrinks to zero as the rotation speed approaches the radial trap frequency, in contrast to the 
finite difference in the non-rotating systems. 
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Strongly interacting two-component Fermi gases pro- 
vide a unique testing ground for the theories of ex- 
otic systems in nature. In atomic Fermi gases, tunable 
strong interactions are produced using Feshbach reso- 
nance P,|3,I3 ■ By sweeping the magnetic field in the Fes- 
hbach resonance experiments, magnitude and nature of 
the two-body interaction strength changes from repulsive 
to attractive. Across the resonance the s-wave scattering 
length a goes from large positive to large negative values. 
The fermionic system becomes molecular Bose-Einstein 
condensates (BEC) for strong repulsive interaction and 
transform into the Bardeen-Cooper-Schrieffer (BCS) su- 
perfluid when the interaction is attractive. The first 
observations of BEC of molecules consisting of loosely 
bound fermionic atoms 0, [E jol initiated a series of ex- 
plorations 0, S H ES Il3| of the crossover between 
BEC and BCS superfluid. Near the resonance, the zero 
energy s-wave scattering length a exceeds the interparti- 
cle spacing and the interparticle interactions are unitarity 
limited and universal. 

One of the hallmarks of superfluidity, be it bosonic 
or fermionic, is the presence of quantized vortices. Re- 
cently, large number of vortices in a Fermi gas along 
the BEC-BCS crossover are indeed observed jl^] and it 
proofs that this strongly interacting system is really su- 
perfluid. Equilibrium properties of vortices in a rotating 
Fermi superfluid described by BCS theory have alr eady 
been the subject of several studies [3, [ll [H, [13, Ell- 
There have been known several outstanding vortex prob- 
lems associated with the fermionic excitations around a 
vortex core in the quantum limit of a superconductor 
under magnetic field |l9j . 

The successful measurements of the breathing mode 
frequencies in a non-rotating Fermi superfluid |20l I2H 
and the availability of the vortex lattices in a rotating 
Fermi superfluid |l3 have stimulated us to study the- 
oretically the collective excitation spectrum of a rotat- 
ing Fermi superfluid containing large number of vortices 
along the BEC-BCS crossover. So far, both theoreti- 
cal I23, lU 01 and experimental [2^ results have been 
obtained for low-energy modes of a rotating, weakly in- 
teracting condensates with large number of vortices. As 
in the case of a rotating bosonic clouds, the frequencies 



of the collective modes of a Fermi superfluid along the 
crossover in presence of the vortex lattices can be mea- 
sured experimentally with high accuracy. 

In this work, we calculate all low energy frequencies of 
the transverse excitations of a strongly interacting Fermi 
superfluid contains vortex lattices along the BEC-BCS 
crossover. In particular, we find that the Fermi super- 
fluid in the BEC regime enters into the lowest Landau 
level at fast rotation. However, the superfluid in the 
unitarity and the BCS regimes occupy many low-lying 
Landau levels even at fast rotation. At fast rotation, 
the radial breathing mode frequency becomes close to 
2u}r {oJr is the radial trap frequency) in all the regimes 
along the crossover, in contrast to the non-monotonous 
behavior of the breathin g m ode along the crossover in 
the non- rotating systems 0| . 

We consider a large number of vortices formed in a ro- 
tating Fermi superfluid confined by a cigar shaped har- 
monic trap potential T4o(?', ^) = {M/2){uj'lr'^+u!lz'^). We 
shall use the concept of diffused vorticity which has been 
successfully used to study the low energy excitations and 
the dynamics of the weakly interacting atomic BEC in 
presence of large number of vortices 22, 23^ J4J- The 
diffused vorticity, V x v = 211, of a vortex lattice is an 
effective coarse-grain averaged description that simplifies 
the formalism by smoothing out the effect of the individ- 
ual vortices. The concept of diffused vorticity is sufficient 
to describe the dynamics at macroscopic distances, larger 
than the average distance between vortices. 

We assume that the system behaves hydrodynamically 
throughout all the regime. If the system is BCS super- 
fluid, then as long as the oscillation frequency is below 
the gap frequency (Ag/?i) needed to break up a Cooper 
pair, this condition is expected to be fulfilled. Within the 
coarse-grain hydrodynamic description, the equations of 
motion for the density and the velocity field which are 
written by the following continuity and Euler equations, 
respectively. 
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~V • [n(r)v]. 



(1) 



2 



and 



2Mv X J7 + Afv X V X V, 



(2) 



where is the density-dependent chemical poten- 

tial Also, V = {h/2M)'V9 is the superfluid veloc- 
ity of a particle of mass M. The terms 2Mv x ft 
and V[(l/2)Mr2^r^] are the Coriolis and the centrifu- 
gal forces, respectively. The equation of state enters 
through the density-dependent chemical potential. We 
assume the power-law form of the equation of state 
as fi{n) = Cri^ . At equilibrium, the density profile 
takes the form no{r) — (/i/C)^/'^(l — f^)'^/^, where 
fi = [M/2){lo1 - ^'^)R\ is the chemical potential, i?^ 
is the radial size and f = r/R±. We are interested to 
study the radial excitations, therefore, we neglect the 
weak harmonic confinement along the z-direction. 

We linearize the density, velocity field and the equa- 
tion of state fJ.{n) around their equilibrium values as 
n = no{r) + Sn, v ^ vq + 6v and /j,{n) — /i(no) -I- 
(9/i/9rt)|„=„p(5n. In the co-rotating frame, the equilib- 
rium velocity field Vq — 0. Thus we obtain the equations 
of motion for the density and the velocity fluctuations in 
the rotating frame as 



d6n _ , / X c- 1 
— = -Vr ■ [no(r)(5v]. 



(3) 



and 



= -Vr[^^\n=noSn] + 2MSv X n. (4) 

We shall follow Ref. to determine the relationship 
between the density and the phase fluctuations in pres- 
ence of the vortex lattices. We are looking for the normal 
mode solutions of the density and phase fluctuations as 
Sn (X e™"^e-*"* and SO oc e™*e-"^*, where m is the az- 
imuthal quantum number, (f> is the polar angle and lu is 
the excitation frequency in the rotating frame. Equation 
Q can be written in terms of the phase fluctuations after 
separating into radial and angular components as 



ihuj- 



dS0 
dr 



dr dn 



^Jn] - 2ihn—6e, (5) 

r 



and 



dii{n) ^ 959 

iriLodO = 2 — - — \n=nnOn — 2inu —. (6) 

on m or 

Integrating Eq. 10 with respect to r, we then get. 



ihojSe = 2 



dy.{n) 



Equating equations © and {7)), one gets 



r or 



6n - 2imhn / dr-SO. (7) 



(8) 



Differentiating the above equation on both sides with re- 
spect to the radial coordinate r and it gives the Cauchy 
equation: 



(9) 



The general solution of the above equation is given by 
69 = c^r'^ + c~r~™' , where c"*" and c~ are arbitrary con- 
stants. The velocity fluctuation along the radial direction 
can be obtained as 5vr = {fi/2M){m/r){c'^r™ — c~r~"^). 
The radial velocity can not diverge as r — > 0. Therefore, 
c~ — for m > and c+ = for m < 0. This gives us 

dS9 ^^^^ 



m , 
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Substituting Eq. lfTU|l into Eq. we obtain the fol- 
lowing relationship between the density and the phase 
fluctuations: 

i(hLu±2hn)S9 ^2^^$^L=„Jn. (11) 
on 

Taking gradient of Eq. Hll|l and then substituting it 
into Eq. Q, we then get, 

= ^(l-f^2)V^.[(l-f2)l/7V_(l_ ^2)1-1/7 



a;ct((I'a ± 2fl)Sn{r), 



(12) 



where uj — uj/ujr and fl — n/tUr- Here, a : {nr,m) is 
a set of two quantum numbers, where n^ is the radial 
quantum number. 

When = 0, the solutions of Eq. H12|l are known 
analytically |23|. For = 0, the eigenspectrum is given 
by 

Col^\m\ + 2nr[-f{nr + \m\) + l], (13) 
and the corresponding eigenfunction is given by 

(5n„(f,0) oc (1 -f2)l/7-l^:|m|pa/7-l>i|)(2f2 _ ly^'l' ^ 

(14) 

where Pn°'''''\x) is a Jacobi polynomial of order n. As- 
suming the net effect of a coarse grain averaged vortex 
lattice on the wave functions to be a modification of the 
effective trapping potential. Therefore, we shall use the 
above wave function as our ansatz. With this ansatz for 
Sn{f, (p), the expectation value of Eq. 112|) is given by 

Wa('I'a±2f2) = [\m\ + 2nr-finr + \m\ + l/-f)]{l-h'^). (15) 

The above quadratic equation can be solved to get the 
following eigenfrequencies in the rotating frame: 



'[\m\ + 2nr-f{nr + |m| + 1/7)](1 - n^) + 

(16) 

The eigen frequencies in the laboratory frame are ob- 
tained by using the simple relation uj-^ot — ^lah ~ tt^^ and 
it is given by 



uJa = \/\\m\ + 2nrj{nr + \m\ + l/"f)]{l-n'^) + 

± (|m| - (17) 
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The frequencies of the m = modes are identical in both 
the rotating and the laboratory frames and the energy 
degeneracy is restored for the m = ±1 modes in the 
laboratory frame. The dipole mode gives wo.ii — 1 for 
all which satisfies the generalized Kohn's theorem. 

When 7 = 1, it reproduces the known result for weakly 
interacting rotating BEC with a large number of vortices 
|24j . The eigenfrequencies given in Eq. (inil with 7=1 
successfully describes the low energy transverse excita- 
tions of a rotating atomic BEC containing large number 
of vortices In particular, the first (wi p) and second 
(^2,o) radial breathing modes are in good agreement with 
the experimental as well as numerical results of weakly in- 
teracting BEC 25, 28]. Moreover, the quadrupole modes 
(^o,±2) are in good agreement with the numerical results 
obtained by solving the Bogoliubov-de Gennes equation 

At zero temperature, the energy per particle of a dilute 
Fermi system can be written as e = {3/5)Epe{y), where 
Ep — h'^kp/2M is the free particle Fermi energy and 
e{y) is a function of the interaction parameter y = 1/kpa. 
The Monte Carlo calculations of Refs. 113,1^ find e(j/ 
0) = 0.42 ± 0.01 which is in fair agreement with the 
experiments 0, On the basis of the data of Ref. 

[8ll |. Manini and Salasnich proposed the following 
analytical fitting formula of e(y) for all the regimes along 
the BEC-BCS crossover: 



e(y) = ai - a2 tan 



-ir /3l 

["32/- 



/32 



\y\ 



(18) 



This analytical expression is well fitted with the data of 
Ref. 31] for a wide range of y on both sides of the reso- 
nance. Two-different sets of parameters are considered in 
Ref. ^3] '■ one set in the BCS regime (y < 0) and another 
set in the BEC regime (y > 0). In the BCS Hmit, the 
values of the parameters [33] are ai = 0.42, 02 = 0.3692, 
as = 1.044, f3i = 1.4328 and P2 = 0.5523. In the BEC 
limit, the values of the parameters 1^3 are ai = 0.42, 
aa = 0.2674, ag = 5.04, /3i = 0.1126 and^ = 0.4552. 
The chemical potential /i is given by |32| 

^ = e{n) + n^ = Ep[e{y) -^e' (y)] = EpF{y), (19) 
an 5 

where e'(y) = de{y)/dy and F{y) = e{y) - (y/5)e'(y). 
Due to the effect of the solid body rotation and the har- 
monic trapping potential, the free Fermi energy can be 
written as Ep — fiujr[iN\{l — Q?)]^/^ with A = Wz/wr- 
The radial and the axial sizes can be written, respectively, 
as Ra_ = a,.(24iVA)i/6(l - n"^)'^/^ ./F{^ and i?^ = 
a^{2AN/\^Y/^{l-h^Y/^./F{y), where a, = ^fi/Muo, is 
the harmonic oscillator length with i = r, z. The aspect 
ratio is given by Rz/Rj_ = (1 — fPy^^/X which shows 
that the rapidly rotating Fermi superfluid takes a pan- 
cake form even if the confining trap has a cigar shape. It 
is interesting to note that the aspect ratio is the same for 
all the regimes along the BEC-BCS crossover. 

It is seen from Eq. H17|l that when ^ ujr all trans- 
verse modes with m = in all the regimes along the 
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FIG. 1: Plots of the chemical potential vs coupling parameter 
y of a trapped Fermi gas with TV = 2 x 10*' and A = 23/57. 



BEC-BCS crossover approach to the degenerate value 
2u}r- Such a macroscopic degeneracy is reminiscent of 
the degeneracy of the energy levels of a single particle 
in an uniform magnetic field, leading to the well known 
Landau level structure. However, we should mention 
that our results here hold not for a single particle, but 
for a system of interacting particles along the BEC-BCS 
crossover. Therefore, the strongly interacting, rotating 
Fermi superfluid forms the Landau level structure when 
Ur- In Fig. 1, we plot the chemical potential for 
various values of fi. Fig. 1 shows that the chemical po- 
tential in the BEC side is less than the cyclotron energy 
gap (2hujr) at the fast rotation and the system enters into 
the lowest Landau level. However, at the fast rotation, 
the chemical potential in the unitarity and in the BCS 
regimes are quite large compared to the cyclotron energy 
gap. The fermionic superfluid in the BCS as well as in 
the unitarity regimes occupies many low-lying Landau 
levels even at — > w^- 

One can extract an effective adiabatic index 7 and its 
dependence on the scattering length a by defining the 
logarithmic derivative as (32j 



^ ^ dn 



2 

_ 3 



fe'iy)+yj,e"{y) 



(20) 



From Eq. H17|l . the first radial breathing mode fre- 
quency can be written as 



cji = y^2(7 + i)(i-fi2) + fi2 + 



(21) 



When il = 0, the difference between the breathing 
mode frequencies at the BEC limit and unitarity limit 
is A(^2 = 0) = a;i(7 = 1) - ^1(7 = 2/3) = (2 - 
■y/l0/3)aj,. ~ 0.174258ti'r- Recently, the radial breath- 
ing mode frequency of a non-rotating Fermi superfiuid 
along the crossover has been measured [23, El which is 
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FIG. 2: Plots of the first breathing mode frequency for S7 = 
(solid line), Cl = 0.75 (dashed line), Q = 0.85 (dot-dashed 
line) and Cl = 0.97 (dotted line). 

in fair agreement with the theoretical result cji(f2 = 0) = 
\/2{j + l)^^^- In Fig- 2, we plot the first radial breathing 
mode in all the regimes along the BEC-BCS crossover for 
various values of fi. It is interesting to see from Fig. 2 



that the difference between the breathing mode frequen- 
cies at the BEC-limit and the unitarity limit decreases 
as ft increases. In fact, the breathing mode in all the 
regimes approaches 2ujr when Q ujr- The shrinking of 
the difference of the breathing mode frequency could be 
observed in the future experiment. 

In conclusion, we have obtained an analytic expression 
for the collective excitation spectrum of a Fermi super- 
fluid containing vortices along the BEC-BCS crossover. 
The spectrum reduces to the Landau level structure in 
the fast rotating regime. The superfluid in the BEC 
regime goes into the lowest Landau level when f2 — > w,.. 
However, the superfluid in the unitarity and the BCS 
regimes occupies many low-lying Landau levels at the 
fast rotating regime. The difference between the breath- 
ing mode frequencies at the unitarity and the BEC limits 
shrinks to zero in the fast rotating regime. 
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